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I. INTRODUCTION 


In the last few years, there has been considerable growth in the subject of modelling 
of materials making use of potent semiempirical methods. These techniques vary in com- 
plexity, physical foundation, numerical simplicity and range of applications but overall, they 
facilitated the development of computational materials science to the level in which it is 
found today. However, several limitations still exist, and much of the recent work done in 
this field is directed toward overcoming these difficulties. This problem is more noticeable 
in the case of alloys: as opposed to the problem of atomic structure of monatomic systems 
where many techniques, including first-principles approaches, have reached a noticeable de- 
gree of sophistication and accuracy, the problem of alloys is, by its own nature, much harder 
to tackle. 

The wealth of experimental studies of surface relaxation on pure metallic surfaces is 
not matched for alloys [1] However, in spite of the small number of experimental studies 
[2-4] there seems to be slow but sure progress in the field, as the available theoretical tools 
for modelling become more accurate. Among them, semiempirical techniques have enjoyed 
widespread use, helped by the ever-increasing computational capability. In this paper, we 
will concentrate on one of these approaches, the BFS method [5], which has shown great 
promise in becoming a viable alternative to deal with the various problems associated with 
alloys, in that it provides an extremelly simple formalism, a sound physical foundation 
based on perturbation theory via equivalent crystal theory [6], and a convenient numerical 
efficiency making it a good candidate for computer simulations of materials properties. 

Since its inception two years ago, BFS has been applied to a variety of problems, starting 
with the basic analysis of bulk properties of solid solutions of fee [5] and bcc [7] binary alloys 
(heat of formation [5], lattice parameter [8], etc.) and more specific applications like the 
energetics of bimetallic tip-sample interactions in an atomic force microscope [9] as well 
as Monte Carlo simulations of the temperature dependence of surface segregation profiles 
in Cu-Ni alloys [10]. Other applications include surface structure [11] and a diagramatic 
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analysis of ordered alloy clusters for the determination of the ground state structure of a 
given binary alloy [12]. An additional advantage of BFS is that it allows for deriving simple, 
approximate expressions which describe the trends in segregation as well as elucidating the 
driving mechanisms for these phenomena [13]. Also, as a consequence of the ideas underlying 
the foundation of BFS, simple expressions for predicting the composition dependence of bulk 
alloy properties based solely on pure component properties have been recently derived (‘BF 
rule’) [14] providing an alternative to the commonly used Vegard’s law [15]. 

In this paper, we concentrate on the application of BFS to the study of surface structure 
of metallic alloys. After a detailed description of the method, we present a simple appli- 
cation to the calculation of surface energies, followed by an extensive set of results on the 
characteristics of surfaces of ordered alloys. These results include an analysis of the surface 
energies as a function of crystal face, composition of the surface and type of structure, as 
well as the prediction of the rippling of the surface due to the different relaxations for atoms 
of different species in mixed-composition surfaces. 
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II. THE BFS METHOD 


A. Basic ideas 

We will introduce the BFS formalism [5,7-12,17,13] with reference to the following ex- 
ample. Consider two pure single crystals: one of atomic species A (lattice parameter a e ) 
and one of atomic species B (lattice parameter af ). This will be the initial state. The 
final state will be a certain alloy A - B with lattice parameter a x . The ‘ideal’ process of 
alloy formation is shown in Fig. 1 . Let us focus on one of the atoms in the A crystal. 
Fig. 2 represents the transformation undergone by this atom: there is a change in geometry 
(the lattice parameter changed from af to a x ) and a change in composition (some of the 
neighbors are now B atoms, denoted by dots). In BFS, we approximate this ‘ideal process 
of alloy formation by a sequence of two independent transformations, as shown in Fig. 3. 

In the first transformation, the identity of the atoms is conserved. The atom in question 
(denoted in figs. 2 and 3 by ®) sees its environment changed only in terms of the relative 
distances of the atoms surrounding it. This is a defect that can be straightforwardly treated 
with any method for monatomic crystals. In particular, because of its effectiveness for 
dealing with this kind of defects, equivalent crystal theory presents itself as a good candidate 
for such calculation. In this transformation, the atom in question suffers a change in energy 
e s , which we will call strain energy, because it is related only to lattice deformations. 

In the second transformation (fig. 3.b), the geometry of the equilibrium crystal is con- 
served. The atom in question sees its environment changed only in terms of the identity of 
its neighbors. For evaluating this change in energy, which we will call chemical energy, we 
assume that the neighbors are located in equilibrium lattice sites. A change in coordination 
in the alloy (for example, atoms close to a vacancy or to a surface) would introduce struc- 
tural information in the chemical energy which, as mentioned above, should be accounted 
for only by the strain energy term. Therefore, we need to reference the chemical energy so 
computed to a similar situation where the neighbors of the atom in question are forced to 
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have the same chemical identity. This way, the chemical energy will only carry information 
on the chemical interaction between atoms of different species, regardless of the geometrical 
distribution or coordination. 


B. Formalism 

Adding the strain and chemical energy contributions for a given atom yields the net 
contribution of this atom to the energy necessary to assemble the alloy A — B from its pure 
components, which is precisely the heat of formation. Therefore, the contribution of atom i 
of atomic species X (X = A, B , ...) to the heat of formation of the alloy is 

e,- = ef -f (ef — ef° )• (1) 

where ef° denotes the reference to the chemical energy. The total heat of formation , AH, 
is just the sum of the individual contributions of each atom in the alloy. In writing eq.(l) 
we assume that the two transformations provide a good representation of the process of 
alloy formation. However, due to the nature of these transformations (i.e., in the calculation 
of the structural energy we ‘freeze’ the chemical composition and in the calculation of the 
chemical energy we ‘freeze’ the atomic locations), the chemical energy is a constant, solely 
dependent on the chemical composition of the alloy, which is clearly unrealistic for some 
situations. One would therefore expect eq.(l) to be a good representation of the alloy for 
those cases where either there is a small lattice mismatch between the pure components or 
when the defect represents a small departure from equilibrium. For those cases where such 
conditions are not met, it is necessary to ‘recouple’ the two independent processes described 
above, weighing the chemical contribution to the heat of formation accordingly. In order to 
accomplish for this, we then introduce a factor which accounts for the asymptotic behavior 
of the enthalpies of formation of alloys. To be defined later, this ‘glue’ factor, <7i, links 
the strain and chemical contributions not only providing a better description of the alloy 
formation process by recoupling the strain and chemical contributions, but also by giving 
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the correct behavior for large interatomic distances. For alloys of elements with small lattice 
mismatch, this term is of little relevance in the region of interest, when the typical distances 
in the alloy are comparable to the distances in the equilibrium crystals. We then rewrite 
eq.(l) as 

£. = ef +9i(ef -ef 0 )- ( 2 ) 


The strain energy, being that it arises from a single-component system can be written as 

ef = E$F*(af*) (3) 

where, based on the assumption that the universal binding energy relation of Rose et al. [16] 
contains all the relevant information concerning a single-component system, 

F*(a*) = 1 — (1 + a*)e _ °*, (4) 


and where is the cohesive energy of a pure crystal of species t. The scaled lattice 
parameter af * is given by 


where af is the lattice parameter of the equivadent crystal associated with the defect, a* and 
li are the equilibrium lattice parameter and scaling length [16] of a pure crystal of species i 
respectively, and q 3 = 3/16* for fee crystals. Although ef can be obtained via ECT [6], it 
does not have to be necessarily so, and af* can be obtained as a solution of eq.(3) if £,• is 
computed by any other technique. Either way, af* can be readily obtained, with which we 

define the coupling term </,- as 


9i 



( 6 ) 


As in previous efforts, we choose ECT [6] to perform strain energy calculations, the choice 
being guided by the simplicity and reliability of this technique. Using ECT for computing 
ef introduces the added advantage that af (and thus af*) is directly obtained by solving the 
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ECT equation for the defect crystal, as shown below. Within the framework of ECT, af is 
interpreted as the lattice parameter of an ideal, perfect crystal (i.e., the equivalent crystal) 
where the energy per atom is the same as the energy of atom t in the actual, defect crystal. 
In general, the ECT equation for computing the strain energy reads 

NR\e~ aRl + M J RSe- (0,+ * ) * 3 = £ r?e-< 0,+s(r ' ))r > (7) 

3 

(see ref. [6] for details) where the quantities p,a, X and the screening function S{r) are 
defined in ref. [6]. The sum on the r.h.s. of Eq. (7) runs over all neighbors of atom i at a 
distance Tj. Eq. (7) is then solved for the lattice parameter of the equivalent crystal af. R\ 
and R? are the corresponding nearest- and next-nearest-neighbor distances in the equivalent 

crystal. The strain energy is then computed with Eq.(3). 

Rigorously, the computation of the strain energy includes four terms (see ref. [6]). In 
this work, we neglect the three- and four-body terms dealing with the bond angle and 
face-diagonal anisotropies and retain only the two-body term that accounts for bond-length 
anisotropies, which we expect to be relevant for atoms in the top (surface) layers. The higher 
order terms would be proportional to the small local fluctuations of the atomic positions 
around the equilibrium lattice sites. We expect that the leading term, Eq. (3), will properly 
account for these small distortions. 

The chemical contribution ef is obtained by an ECT-like calculation. As opposed to the 
strain energy term, the surrounding atoms retain their chemical identity , but are forced to 
be in equilibrium lattice sites. If N ik (M ik ) denotes the number of nearest(next)-neighbors 
of species k of the atom in question (of species t) then the ECT equation to be solved for 
the equivalent lattice parameter af is 

NRFe- aiRl + MR^e~ (ai+ ^ )R2 = £ N ik rfe~ aiiri + £ M ik ife~ {a,k+ ^ )T2 (8) 

Jfc k 

where N(M) is the number of nearest(next)-neighbors in the equivalent crystal of species 
i and Ri(R?) is the nearest(next)-neighbor distance in the equivalent crystal of lattice pa- 
rameter af. ri and r%, are the equilibrium nearest- and next-nearest- neighbor distances in 
an equilibrium crystal of species t, respectively. The chemical energy is then computed with 
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ef = 7 B- 0 r(*f-) 


(9) 


and 

«g = (10) 

where 7(70) = +1 if af’(a%*) > 0 and 7(70) = -1 otherwise, and af* = q(af - a‘ e )/4 
(af°* = q(af° — a‘)/4). The scaled lattice parameter af* is obtained from Eq.(8) with the 
parameters aik listed in ref. [17], and af°* is computed by solving Eq.(8) but with a,-* = o,-. 

C. The BFS parameters A 

Using eq.(8) for the calculation of the chemical energy might suggest that this is a 
rigorous ECT calculation. It is not. In interpreting the change of composition as a ‘defect’, 
we are just adapting the basic concept underlying equivalent crystal theory to this case. 
In the single-crystal ECT, where all the atoms are of the same atomic species, we apply 
perturbation theory in order to find the energy of the defect crystal. The perturbation is 
basically due to the difference in potentials between the defect solid and the ground-state 
crystal. As described in the original formulation of the method [6], it is reasonable to 
parameterize the first-order contributions to the perturbation expansion as 

A E oc R p e~ aR (11) 

where p = 2n — 2 (where n is the atom principal quantum number) and a is a parameter 
that will primarily reflect the structure of the electron density in the overlap region. In 
single-crystal ECT the parameter a is determined for metals so that the energy to form a 
rigid (unrelaxed) vacancy is equal to the experimental value. 

To a good approximation, these concepts should remain valid in the case of alloys, and 
by using eq.(8) we are adopting the same functional form used in ECT to describe the 
perturbation due to dissimilar atomic species. In order to deal with arbitrary defects and 
structures in future applications, as well as with multicomponent systems, it is convenient 
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to ‘localize’ this effect and assume that the global property parameterized by a (i.e., the 
tails of the overlapping electron densities) can be separated into pairs of interacting atoms. 
In this approximation, the electron density in the region between two atoms of the same 
species would not be affected by the presence of neighboring atoms of different species. The 
perturbation would then be localized in the region between two dissimilar atoms. This 
assumption justifies the definition of the parameter a,* as 

or,* = a, - + A*,- (12) 

where a, is the usual value of a for the pure element t and A*,- is a correction introduced 
by the presence of a neighbor of species k. Obviously, A*,- = 0 if t = k. The ‘perturbation’ 
parameters A*, and A,* are the only new parameters introduced in BFS: all other parameters 
are those corresponding to the pure components of the alloy. Generally, these two parameters 
are determined by requiring BFS to reproduce the experimental values of the heat of solution 
in the dilute limit, E BA (the heat of solution of an impurity B in a host A), given by 

(13) 

x -+0 

where A H is the heat of formation of the compound A\- X B X , x being the concentration of B 
atoms. These parameters have been computed for a variety of fee and bcc metallic alloys (see 
ref. [17] for a discussion) and work is under way for hep-based structures. For completeness, 
yyg Include both the ECT and BFS parameters for some fee and bcc elements. Table 1 lists 
the ECT parameters p, l, a, A, the cohesive energy Ec and the lattice parameter, while 
Table 2 displays the values of the BFS parameters Aab and Aba for some of the alloys of 
these elements, as well as the experimental values of the heat of solution (eq.(13) in the dilute 
limit used for fitting the BFS parameters. For the case of Al-Ni, we used the theoretical 
values computed by Sanchez and Carlsson [18]. 
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III. SURFACE ENERGIES OF ORDERED STRUCTURES 


A. Ordered structures of fcc-based binary alloys 

In this section we will apply the BFS formalism to the simple problem of computing 
the surface energy of rigid ordered alloy structures. By rigid we mean that, for the sake 
of simplicity, we will not allow for individual or collective atomic displacements due to the 
presence of the surface. However, later on we will show results involving relaxations and 
compare them with available experimental data. 

Several alloys form simple ordered structures for low temperatures [19]. Obviously, the 
higher the symmetry in the patterns that characterize a given structure, the simpler the 
calculation is. Fortunately, many alloys are found to have, for specific concentrations, some 
of the simplest possible structures. Therefore, we will concentrate on these structures, 
shown in fig. 4, for three concentrations: the Ll 2 A 3 B or AB 3 structures and the Ll 0 and 
Lli structures at 50 % concentration, both being fcc-based. Within the cluster expansion 
method [21] these structures, together with the pure fee A and B crystals, form a set of 
fundamental structures for fee alloys as these represent the possible structures that can 
be formed assuming nearest-neighbors interactions only, where a tetrahedral cluster is the 
unit cell [22]. Therefore, the clusters A m B 4 - m (m = 1,...,4) axe the building blocks of the 
corresponding fee alloy structures Ai, Ll 2 and LIq. 

B. A simple example 

We will consider a simple example: the calculation of the surface energy of a given termi- 
nation of a Ll 2 structure of an fee binary alloy. There are two possible terminations for the 
(100) Ll 2 structure of a certain A 3 B alloy: a mixed-composition (1:1 A:B) plane alternating 
with a pure A (1:0 A:B) plane, giving an overall stoichiometry (3:1 A:B), and a pure A 
plane alternating with mixed-composition planes. These two possible bulk truncations are 
also possible for the (110) surface, whereas the (111) truncation is always stoichiometric (3:1 
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A:B). We will concentrate on the (1:1 A:B) Ll 2 (100) surface and label the planes parallel 
to the surface with the index j = 1 , b where j = 1 corresponds to the surface plane and b 
labels a certain bulk plane, below which no surface effects are to be considered. By forcing 
the semi-infinite slab to be rigid and limiting BFS to deal with second-neighbor interactions, 
the calculation is greatly simplified, as only two planes (j = 1, 2) contribute to the surface 
energy. Because of the high symmetry of this structure, only one atom of each species needs 
to be considered for the j = 1 plane, whereas only one A atom contributes from the plane 
j = 2. Atoms in layers j = 3, ..., b are inert, in that they will not contribute to the surface 
energy. We will assume that the lattice parameter a is the one previously determined for 
a bulk ordered alloy A Z B. Let X, ( X = A, B; j = 1,2 ) denote a non-equivalent atom of 
species X in layer j. We need to compute the contributions c Xj from Ai, B\ and A 2 so that 
the surface energy will be 

<T = + GB\ + 2<Tjta (14) 

where 

*Xj = £ x, + 9Xj( £ Xj ~ £< x } ) 

Using eq.(7), the strain energy contribution is computed from the solution of the following 
equations: 

+6;^e■ < °* + £ )R> = 8r^*e-“■•'‘+5r5• , e' < "•‘ + £ ,, ’ (X,=A,) (16) 

= (X, = fl,) (17) 

and 

12 R[ A e~ aARl + 6B% A e~ (aA+ '^ )R * = 12^e-^ ri + 5r$ >l e' (oA+ £ )r8 ( Xj = A 2 ), (18) 

where r*i = \/2a/2 and r 2 = a. The equivalent lattice parameters a^. are obtained from the 
solutions to eqs.(16)-(18) as a s Xj = \/2f?i(j) = Mi) the strain ener gy contributions 
e Xj are computed using eq(4). 
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Using eq.(8), the chemical energy contribution is computed from the solution of the 
following equations: 

12R\ A e~° ARl + 6R% A e~ { ° A+ ^ )R2 = 4<% A e~ aAdl +4(iJ i, e‘ MAw)Jl +5<J? 2 A e~ ( ° A+ ^ )d2 

( Xj = A,) (19) 

12 R\ B e~ aBRl + 6R% B e~ {QB+ ^ )R2 = 4 s[ B e~ aBl1 + 4s^ B e- (ai,+A>u,),, + 5s$ B e" (afl+ ^ ),a 

(Xj = B l ) (20) 

and 

l2R{ A e~ aARl +6R^ A e~ {QA+ ^ )R2 =Sd p 1 A e~ aATl +i<R 1 A e~^ A+ABA)dl +br^ A e~ (aA+ ^ )T2 

(Xj = A 2 ), (21) 

where d x = v/2a?/2, r 2 = a?, s a = \/2a? / 2 and s 2 = af . The equivalent lattice parameters 
a,x are obtained from the solutions to eqs.(19)-(21) as af^. = y/2 Ri(j) = Ri(j)- Remem- 
bering that lx, = ±1 whether a£* = q(a% - a*)!h is positive or negative, the first term 
in the chemical contribution, is then 

t% = ( 4 ‘) < 22 > 


Finally, the reference to the chemical energy, e*° is computed from the solution of the 
following equations: 


12 Rl A e~° AR ' + 6RZ A e~ iaA+ ^ )R2 = + §£ 2 A e~ {aA + 

£ 

II 

>-• 

(23) 

12R{ B e- aBR ' +6R^ B e~ (aB+ ^ )R2 = 8 + 5s5 B e _(aB+ ^ ) ' J 

(Xj = Bi) 

(24) 

and 



12i2f x e- a ^ +6i^ >l e _( “ yl+ ^ )fi2 = \2dF x A e~ aAdl ->rb^ A e~ {aA ^ )d2 

(Xj = As), 

(25) 

The equivalent lattice parameters a£° are obtained from the solutions 

to eqs.(23)-(25) as 


a Co _ yfeR^j) = i? 2 (j) and the reference chemical energy contributions are computed 
using eq(10). 
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IV. APPLICATION TO SURFACE STRUCTURE 


The determination of the surface structure of alloys has recently been the topic of theoret- 
ical and experimental work. Much effort has been devoted to the study of surface relaxation 
in metals and, to a lesser extent, alloys. Several recent experiments have provided insights in 
the phenomena of surface relaxation and composition, in the case of alloys, and correspond- 
ingly a number of theoretical studies have shown good general agreement with experimental 
results. 

However, there is still a great deal of uncertainty in certain areas, due to limitations 
inherent in experimental techniques and also to the lack of alternative studies to verify 
previous results. 

The first experiment, in 1984, that provided detailed information on the atomic positions 
of surface atoms in a truncated ordered alloy is the low-energy electron diffraction (LEED) 
intensity analysis of Davis and Noonan of a NiAl(llO) surface [2]. They found strong 
evidence for a rippled surface, where the A1 sites of the top layer (in the mixed-composition 
truncation) are displaced above the Ni sites by approximately 0.22 A. This result was quickly 
followed by the calculation performed by Chen, Voter and Srolovitz using the embedded atom 
method (EAM) [23], which confirmed the main features found in the experiment. EAM 
was later used to investigate similar phenomena in other ordered alloys: Foiles and Daw 
presented a complete study of Ni 3 Al (Ll 2 structure) [24], followed by Foiles work on ordered 
surface phases of Au on Cu [25], and Lundberg’s extensive study of surface segregation and 
relaxation of Pt-Ni alloys [26]. At the same time, new experimental LEED results on Ni 3 Al 
were reported by Sondericker and coworkers [3], finding a similar rippled structure in Ni 3 Al 
(100) faces. Finally, a low-energy ion scattering spectroscopy (LEISS) experiment by Wang 
and coworkers provided similar data for the Cu 3 Au system, a classic ordering alloy [4]. Their 
work followed the LEISS results concerning the surface composition of the top atomic layers 
[27]. This system was also the subject of a very recent study by Wallace and Ackland using 
a molecular statics algorithm with Finnis-Sinclair (FS) many-body potentials [28]. 
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In this section we will show new results concerning the above mentioned systems as 
well as predictions for other alloys, using BFS. Due to its computational efficiency, this 
application of BFS to the surface structure of ordered alloys extends to a number of systems 
for which there are no other theoretical or experimental studies to date. 

Before proceeding to the calculation of multilayer relaxation in alloys, we will discuss 
some features of theoretical calculations of these quantities. Ref. [1] provides a reasonably 
large sample of both experimental and theoretical results for changes in interlayer spacing 
in pure crystals. In all cases, the theoretical techniques used rely either on the use of 
input data (generally experimentally determined) or on certain approximations for some of 
the variables of relevance. Necessarily, the results will depend on such choices. Multilayer 
relaxations involve at best very small changes in position, and correspondingly, comparable 
changes in surface energy, whose minimization is the criterion used to determine the final 
interlayer spacings. Thus, the search for a minimum of the surface energy, as accurate as the 
minimization technique might be, will be strongly influenced by the two factors indicated 
above: the approximations used and the shallowness of the minimum in the surface energy 
surface resulting from small changes in the input parameters. As a consequence, to quote 
just one value for each of the changes in interlayer spacings as is ordinarily done, might 
not reflect the ambiguities in these calculations. In this paper we adopt a different path: 
to each theoretical prediction, we will attach an estimate of the possible errors due to 
any of the reasons mentioned above. Although there is no certain way to determine such 
errors (after all, the predictions are, within their own framework, exact), we will see that 
changes on the order of 1 % in the surface energy can generate quite interesting variations 
in the relaxation schemes predicted. In particular, within the framework of ECT, such small 
changes in the surface energy can be easily obtained by changing any of the input parameters 
(lattice constant, cohesive energy, bulk modulus) by a similar amount, well below the usual 
experimental errors in the determination of such quantities. 

To illustrate this issue, we will focus our attention on the surface structure of some 
fee pure metals (Al, Au, Cu and Ni). As can be seen in Tables 2-11 of ref. [1], previous 
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theoretical and experimental studies show a wide spread in the predictions of the changes 
in interlayer spacings for the (100) and (110) surfaces. Even results obtained within the 
same theoretical te chni que (EAM, ECT) do not agree with each other (due to different 
fitting procedures of the embedding function in the case of EAM and different input data 
in both cases). Although there is general qualitative agreement, regarding the contraction 
or expansion pattern found for successive layers, in some cases the theoretical values show 
poor agreement with experimental results (see, for example, A1 (100)). The ECT results 
(from refs. [1] and [6]) also highlight this inconsistency. The difference between the values 
obtained in this work and those from previous applications of ECT is easily traceable to 
slightly different values of some of the input parameters. 

As mentioned above, in order to account for these and other ambiguities in the calcula- 
tion, we investigated the change in predicted relaxations due to small changes in the rigid 
surface energy. We thus defined t error bars* in such way that all the intermediate values so 
obtained predict variations in surface energies within that tolerance. Needless to say, this 
range of values does not include all the possible sets (Adu, Adn) that correspond to surface 
energies within the allowed values. It is interesting to note, however, that in most cases, all 
the experimental as well as theoretical predictions fall within the range of values obtained 
in this fashion. 

It should be noted that when comparing our theoretical predictions with available experi- 
mental results, the error bars quoted in each case are not rigourously comparable. However, 
we choose to do so with the only purpose of giving a complete description of the results 
obtainable with ECT and BFS, once uncertainties in the input parameters are taken into 
' account. To illustrate this point, we first discuss the surface energies and multilayer re- 
laxations of the unreconstructed low-index surfaces of pure Al, Ni, Cu and Au crystals. In 
Table 1 we display the ECT predictions for the surface energies and compare the results with 
typical experimental values for polycrystalline samples [29-31]. The agreement is excellent 
in all cases, with the exception of Cu, that shows a somewhat larger deviation than the other 
elements. We note that experimental values for the surface energies are for polycrystalline 
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surfaces, thus could be strongly dominated by the predominant surface plane. 

In table 2 we compare results for the multilayer relaxations of the first two interlayer 
spacings for those cases for which recent experimental data is available [32-41]. Once again 
the agreement is excellent, as it was shown in previous ECT studies of surface structure [1]. 
The inclusion of the theoretical ‘error bar’, as mentioned above, allows for a better compar- 
ison with experiment as it shows that for most cases, small changes in the input parameters 
of the method suffice to account for the whole range of possible experimental results. The 
exceptions are Al(100) and Al(lll), where the outward relaxation of the surface layer has 
been attributed to an electron promotion effect [41]. Semiempirical methods (ECT, EAM. 
etc.), unless specifically designed to do so, do not generally allow for such fine electronic 
structure effects, thus it is not surprising that our results for Adi 2 in these cases predict 
surface layer contractions, even when the ‘error bar’ is taken into account. For completeness 
we also include results for the surface relaxation when only that plane is allowed to relax, 
in order to single out the influence of subsequent interlayer spacing changes on the surface 
plane. Again, the agreement with available experimental data is very good in all cases. 

As mentioned above, there are few theoretical or experimental studies of ordered alloy 
surfaces [2-4,23-28]. First, we discuss two cases (NhAl and C113AU, in the LI2 structure) 
which have been the subject of recent studies [3,4,23,24]. We follow this discussion with a 
complete presentation of the corresponding results for a larger number of ordered structures 
as well as different binary alloys of fee elements, for which no theoretical or experimental 
data exists. 

a. Cu 3 Au: Table 5 displays the results for the unrelaxed and relaxed surface energies 
(in ergs/cm 2 ) as obtained with BFS and with FS many-body potentials [28]. Both meth- 
ods predict, as expected, lower surface energies for the mixed-composition (100) and (110) 
truncations. This feature has been experimentally proven via a low-energy ion scattering 
study which detected equal parts of Cu and Au in the top layer [27]. ECT and FS results 
also agree on the relative change in surface energy once the top-most layers are allowed to 
relax, in spite of the fact that the FS values are 50 % smaller than the BFS ones. As is 
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also to be expected, the surface energies of (100)1:0 and (110)1:0 faces are comparable to 
the corresponding values for single Cu crystals. The corresponding relaxations are quoted 
in table 4. In order to avoid ambiguities in determining the exact atomic positions from the 
entries in table 4, we present the relaxations as the percentage change in interlayer spacing 
from the unrelaxed case to the one measured from the relaxed position to the unrelaxed 
location of the plane immediately below. We also include the BFS predictions for the pure 
Cu truncations of the (100) and (110) planes. Although it is to be expected that the top 
layer relaxation will change as deeper layers are allowed to relax, any ensuing changes would 
be small, not affecting the conclusions drawn from our results. 

For the CU3AU (100) 1:1 Cu:Au case, the results in table 6 imply a rippling of 0.148±0.025 
A, which amounts to 3.97 % of the lattice parameter determined for this alloy (3.73 A). 
This result compares very well with the 3.77 % rippling (A1 out, Ni in) obtained using FS 
potentials. A similar situation is found for the (110) 1:1 Cu:Au surface, where we find the 
rippling to be 4.2±1.1 % of the lattice parameter, whereas FS potentials predict a rather 
smaller change of 1.9 %. For the mixed-composition (111) 3:1 Cu:Au surface, BFS predicts 
a rippling of 4.6±0.4 % thus agreeing with FS results and experimental evidence that the 
Au atoms are farther out than the neighboring Cu atoms in mixed-composition surfaces. 

b. NiqAI The surface energies of relaxed (100), (110) and (111) surfaces are shown in 
Table 7, where we compare our results with the EAM study of Foiles and Daw [24]. As 
found for the Cu 3 Au case, the mixed-composition truncations always have a lower surface 
energy. The differences between the EAM and BFS predictions are consistent with previous 
calculations for pure metals, where the EAM results are typically 50 % lower than the 
experimental ones. Surface relaxations are indicated in Table 8, using the same format and 
notation of table 4. From these results we extract the following values for the gap between 
Ni and A1 atoms in the mixed-composition (100), (110) and (111) surfaces: 0.12±0.04 A, 
0.09±0.05 A and 0.16±0.03 A, respectively. A similar trend, but with somewhat smaller 
values for the rippling are obtained from EAM [23]: 0.09 A, 0.06 A and 0.07 A. A different 
EAM calculation [25] predicts a 0.06 Aseparation between Ni and A1 atoms in all three 
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surfaces. Recent LEED data [3] show A d u (Ni) = -2.73% and A d 12 (Al) ~ 0 (i.e., A1 is 
displaced outward with respect to Ni) and a rippling of 0.02±0.03 Afor the (100) surface 

c .Other fee binary alloys Table 9-11 display results similar to the ones described previ- 
ously for Ni-Al and Cu-Au alloys, indicating, when appropriate, the distance between the 
atoms of different species in the top layer of each surface. Table 9 includes the results for the 
three possible terminations of the Ll 0 structure (1:1, 0:1 and 1:0) whereas Table 10 shows 
the results for the (100)1:1 surface of the Ll, ordered structure. Finally, Table 11 shows an 
extensive set of results of A Z B and AB 3 alloys in the Ll 2 structure, as shown in fig. 4. 

In general, surface energies vary widely for each system, depending on the type of struc- 
ture and crystal face considered. There is no apparent pattern that describes trends among 
all the different systems discussed in this work. For example, although there is some corre- 
lation between the lattice mismatch and the ripple between the surface atoms in the mixed 
termination surfaces of the Ll 0 alloys, in that the vertical separation between atoms of dif- 
ferent species on the surface increases with the lattice mismatch, there is also an exceptional 
case, like Cu-Pd, with a very large rippling effect. The ordering of the surface energies varies 
from alloy to alloy covering a wide range of values. With very few exceptions, the planar 
relaxations are small, with even smaller uncertainty factors which in turn translates into 
small uncertainties in the rippling of the mixed-termination surfaces. 

Obviously, not all the listed structures exist in nature. We are presenting this survey 
with the only purpose of highlighting the simplicity of alloy surface energy calculation with 
BFS and the predictive power of the method. In this sense, the results presented in this 
work will hopefully motivate the development of new experiments in the determination of 
surface structure. 
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V. CONCLUSIONS 


In this paper we presented an extended discussion on BFS, a semiempirical method 
for alloys, concentrating in aspects of surface structure. The application of BFS to the 
calculation of the surface energy of an ordered alloy highlights the computational simplicity 
of this method, which is an essential ingredient of any algorithm used for the determination 
and simulation of materials properties. The results on multilayer relaxation of ordered 
structures, and their comparison with available experimental values, provide a good example 
of the quality of the results that can be obtained with this method. This, together with the 
numerical efficiency of this technique, makes BFS an appropriate tool for dealing with more 
complex situations, like those found in realistic problems in alloy structure. 
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Element 

Cohesive 

Energy 

Lattice 

Constant 

P 

l 

a 

A 

A1 

3.34 

4.05 

4 

0.336 

2.105 

0.944 

Cu 

3.50 

3.615 

6 

0.272 

2.935 

0.765 

Ag 

2.96 

4.086 

8 

0.269 

3.337 

0.756 

Au 

3.78 

4.078 

10 

0.236 

4.339 

0.663 

Ni 

4.435 

3.524 

6 

0.270 

3.015 

0.759 

Ir 

6.94 

3.84 

10 

0.235 

4.408 

0.661 

Pd 

3.94 

3.89 

8 

0.237 

3.612 

0.666 

Pt 

5.85 

3.92 

10 

0.237 

4.535 

0.666 

Fe 

4.27 

3.57 

6 

0.279 

2.963 

0.784 

W 

8.66 

3.16 

10 

0.274 

4.232 

0.770 

Ta 

8.10 

3.30 

10 

0.325 

3.905 

0.914 

Mo 

6.82 

3.15 

8 

0.262 

3.420 

0.736 

Nb 

7.57 

3.30 

8 

0.341 

3.243 

0.958 

V 

5.31 

3.03 

6 

0.305 

2.726 

0.857 

Cr 

4.10 

2.88 

6 

0.254 

2.889 

0.714 

Fe 

4.29 

2.86 

6 

0.277 

3.124 

0.770 

Li 

1.63 

3.491 

2 

0.589 

1.049 

1.66 

Na 

1.113 

4.225 

4 

0.578 

1.359 

1.62 

K 

0.934 

5.225 

6 

0.694 

1.528 

1.95 

Rb 

0.852 

5.585 

8 

0.651 

1.937 

1.83 

Cs 

0.804 

6.045 

10 

0.757 

2.115 

2.13 


Table 1: Experimental input: Cohesive energy (in eV), lattice parameter (in A). ECT 
parameters: p, l (in A), a (in A -1 ) and A (in A) for several fee elements. 
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A-B 

&AB 

&BA 

IT 1 JO/ 

&AB 

Wa 

Ag-Al 

0.0475 

-0.0499 

0.141 

-0.166 

Ag-Au 

-0.0333 

-0.0227 

-0.161 

-0.186 

Ag-Cu 

-0.0391 

-0.0308 

0.392 

0.250 

Ag-Pd 

-0.0451 

-0.0178 

-0.108 

-0.288 

Al-Au 

-0.0501 

-0.0853 

-1.26 

-0.80 

Al-Cu 

-0.0526 

-0.0626 

-0.35 

-0.20 

Al-Ni 

-0.0657 

-0.0861 

-1.715 

-0.494 

Au-Cu 

-0.0513 

-0.0604 

-0.191 

-0.126 

Au-Ni 

-0.0506 

-0.0622 

0.280 

0.218 

Au-Pd 

-0.0460 

-0.0345 

-0.195 

-0.356 

Cu-Fe 

0.0495 

0.0638 

0.412 

0.349 

Cu-Ni 

-0.0163 

0.0309 

0.100 

0.032 

Cu-Pd 

-0.0431 

-0.0495 

-0.392 

-0.436 

Cu-Pt 

-0.0585 

-0.0441 

-0.299 

-0.532 

Fe-Ni 

-0.0106 

-0.0320 

-0.218 

-0.079 

Fe-Pd 

-0.0229 

-0.0584 

-0.656 

0.177 

Ni-Pd 

-0.0396 

-0.0478 

-0.088 

0.057 

Ni-Pt 

-0.0609 

-0.0537 

-0.330 

-0.282 

Cr-Fe 

0.0465 

0.0285 

0.218 

0.218 

Cr-Mo 

-0.02447 

-0.0090 

-0.102 

0.807 

Cr-V 

-0.0246 

-0.0232 

-0.088 

-0.189 

Fe-V 

0.0998 

-0.07168 

0.215 

0.323 


Table 2: Parameters A ab and &-BA (in A -1 ) and heats of solution in the dilute limit EJ£b 
and (in eV/atom) for several fee binary alloys. The values for Al-Ni were obtained 
from a theoretical calculation by Sanchez and Carlsson [18] . 
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Technique 

Al 

Cu 

Ni 

Exp. [29] 



2270 




2380 




2240 

Exp. [31] 

1169 


2664 

ECT(IOO) 

1203 

2309 

2982 

ECT(llO) 

1284 

2373 

3073 

ECT(lll) 

856 

1767 

2274 


Table 3: Surface energies (in ergs/cm 2 ) of fee Al, Cu and Ni 


Element 

Face 

Experiment 


ECT 

ECT 

(two- layers) 



Adi2 

A<^23 

Ref. 

Adi2 

Adi2 

Ad23 


(100) 

+1-8 


[32] 

-4.68+1.62 

-5.05+1.58 

+3.35+0.80 

Al 

(110) 

-8.5+1.0 

+5.5+1.1 

[33] 

-8.29+2.35 

-9.53+3.58 

+1.90+2.24 


(111) 

+1.7+0.3 

+0.5+0.7 

[34] 

-3.67+1.21 

-3.94+1.19 

+2.75+0.61 


(100) 

-3.2+0.5 


[35] 

-3.53+1.68 

-3.82+1.68 

+2.48+0.85 

Ni 

(110) 

-9.0+1.0 

+3.5+1.5 

[36] 

-6.32+2.44 

-6.55+3.63 

+0.34+2.24 


(111) 

-1.2+1.2 


[37] 

-2.89+1.29 

-3.10+1.25 

+2.12+0.63 


(100) 

-2.1 

+0.45 

[38] 

-3.52+1.74 

-3.81+1.70 

+2.47+0.86 

Cu 

(110) 

-7.5+ 1.5 

+2.5+1.5 

[39] 

-6.31+2.46 

-6.51+3.83 

+0.29+2.44 


(111) 

-0.7±0.5 


[40] 

-2.88+1.30 

-3.10+1.25 

+2.12+0.63 


Table4 : Surface relaxations of Al, Cu and Ni as percentages of the rigid interplanar spacings. 

































mMm 






896 

865 

2119 

1810 


1192 

1171 

2478 

2247 


1051 

1024 

2397 

2337 


1240 

1173 

2873 

2699 

(111)3:1 

882 

863 

1626 

1577 


Table 5 : Unrelaxed and relaxed surface energies (in ergs/cm 2 ) of CU 3 AU. 


Layer 

Atom 

( 100 ) 1:1 

( 100 ) 1:0 

( 110 ) 1:1 


(111)3:1 

H 

Cu 

Au 

-2.12±0.41 

+5.82±0.81 

+5.38±0.47 

-2.45±1.14 

+7.80±1.80 

-3.70±1.02 

+14.35±0.26 

+22.29±0.53 

2 

Cu 

Au 

+5.61±0.37 

+12.87±0.50 

+3.23±1.23 

+10.14±1.37 

EnwgjiM 

MEBm 

+21.20±0.24 

+2.23±1.92 


Table 6 : Surface relaxations of CU 3 AU LI 2 surfaces. See text for definition of the percentage 
change. 
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Table 7 : Surface energies of N13AI (in ergs/cm 2 ). 


Layer 

Atom 

(100)1:1 

(100)1:0 

WB5SSSM 

(110)1:0 


1 

Ni 

-1.33±0.70 

+3.41±0.72 

-4.63±1.65 

-5.03±1.31 

+8.56±0.55 


A1 

+5.25±1.60 

- 

+0.80±2.85 

- 

+16.28±1.40 

2 

Ni 

+4.59±0.63 

+8.82±0.81 

+3.55±2.19 

+5.19±2.40 

+13.69±0.53 


A1 

- 

+3.39±1.46 

- 

+0.92±2.74 

+4.04±1.53 


Table 8 : Surface relaxations of N^Al LI2 surfaces. See text for definition of the percentage 
change. 
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Al-Ni 



Table 9 : Surface energies (in ergs/cm 2 ), planar relaxations (as percentages of the rigid 
interplanar spacing, see text) of the top two layers of several low-index faces of the Llo 
ordered structure. A- indicates the relaxation between layers i and j for an atom of 
species X in layer t. The last column indicates the difference (in A) between the position 
of an atom A and an atom B in the top layer. 
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Au-Pd 


Face 

<T 



^12 

A& 

A A 
^23 

Tv 3 

n 23 

Ripple 

1 1 wins 

1953.39 

-0.37±0.57 

-4.06i0.60 

-0.25i0.50 

+2.76i0.71 

0.07415±0.02355 

' I ffi|| 

1404.01 

-2.84±0.44 



-0.09±0.41 


i 

2668.31 


-1.42i0.80 

+0.32i0.79 



1 1 

1976.47 

-2.96±0.95 



+0.39il.63 


■ ! 

2364.57 


-7.79il.24 

-l.38il.ll 




Nil 

5 d 

I Face 

a 

A* 

n 12 

Ag 

^23 

A& 

Ripple 


1717.93 

-8.55i0.26 

-13.06i0.68 

-7.92i0.30 

+2.44i0.28 

0.08520±0.01768 


2176.18 

-15.98±0.29 



-7.94i0.23 



1255.65 


-9.01i0.19 

-0.71i0.21 




1980.87 

-12.61i0.31 



-5.14il.51 



1902.53 


-19.89±1.16 

-10.47i0.49 




Ni-1 

Pt 

Face 

a 

A& 

Ag 

A23 

A& 

Ripple 

i am b 

2153.11 

-4.63i0.18 

+2.08±0.44 

+5.87i0.19 

-0.87i0.41 

0.12808±0.01181 

' Wm 

3248.61 

-3.75i0.68 



-1.60i0.73 


' i Wro 

2228.15 


-5.60±0.35 

-1.94±0.33 



] I 

2519.11 

-11.94i0.39 



-3.38i0.61 


IHIipb 

3057.54 


-2.66il.13 

-9.56i0.83 




Table 9: (Cont’d) 
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Alloy 

a 

Aj 2 


A 23 


Ripple 

Al-Ni 

1442.33 

+1.97±0.67 

-5.14±0.16 

-0.89±0.54 

+5.18±0.17 

0.13379±0.01566 

Cu-Au 

1205.58 

-2.61±0.15 

+5.25±0.44 

+8.26±0.15 

+0.44±0.41 

0.15463±0.01155 

Cu-Pd 

1816.93 

-2.30±0.39 

-13.64±0.67 

-6.08±0.75 

+8.38±0.41 

0.21487±0.02022 

Cu-Pt 

1964.02 

-3.31±0.24 

+2.81±0.44 

+6.59±0.24 

+0.25±0.45 

0.11756±0.07319 

Ag-Au 

1526.48 

-5.91±0.54 

-4.00±0.35 

+0.69±0.47 

-1.96±0.46 

0.03921±0.01829 

Au-Pd 

1885.63 

-0.47±0.55 

-4.40±0.58 

-0.41±0.49 

+2.89±0.68 

0.07898±0.02290 

Ni-Pd 

1549.94 

-8.06±0.23 

-13.20±0.62 

-8.03±0.27 

+3.04±0.25 

0.09749±0.01626 

Ni-Pt 

1856.99 

-5.02±0.16 

+1.64±0.39 

-5.61±0.16 

-1.18±0.38 

0.12788±0.01051 


Thble IQ : Surface energies a (in ergs/cm 2 ), planar relaxations (as percentages of the rigid 
interplanar spacing, see text) of the top two layers of the (100)1:1 surface of the Llj ordered 
structure. A* represents the relaxation between planes i and j for an atom of species X 
in layer t. The last column indicates the difference in position (in A) between an atom A 
and an atom B in the top layer. 
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CU3AU 


Face 

a 

Aj2 

A 5 

^12 

A 23 

X s 

^23 

Ripple 

n <1*1 lid 

1922.98 

-6.51±0.30 

+1.42±0.66 

+1.67±0.28 


0.14941±0. 01823 

R III 

2205.44 

+2.36±0.33 


+10.47±0.35 

+1.33 ±0.85 


R jm fin 

2453.12 

-3.63±0.85 

+6.58±1.48 

+10.00±0.85 



R mmm 

2649.79 

-8.07±0.61 


+7.56±1.17 

-2.65±0.86 


(111)3:1 

1309.49 

+8.60±0.22 

+ 16.68±0.52 

+15.93 ±0.20 

+2.48±0.87 

0.17583±0.01601 


AU3CU 


Face 

a 

A* 

“12 

VB i 

^12 

A* 

^23 

^23 

Ripple 

Ijjgjgl 

2050.03 

1285.48 

+1.47±0.83 

-0.12±0.26 

-3.89±0.41 

+1.18±0.37 

+0.91±0.48 

+5.03±0.31 

0.10733±0. 02479 

i'll 

1869.36 

1905.90 

-4.02±0.93 

-1.66±0.63 

-6.92±0.64 

+0.44±0.92 

-1.09±1.57 

+8.66±1.06 

0.04096±0. 02221 

(111)3:1 

1286.93 

-2.23±0.15 

-7.73±0.17 

-1.88±0.16 

+3.83±0.17 

0.12722±0.00743 


Cu 3 Pd 


I Face 

O 

A12 

Af 2 

Aj3 

sg 

Ripple 


2341.92 

2259.46 

-5.58±0.51 

-7.17±0.51 

-6.83±0.62 

+0.43±0.36 

-5.89±0.75 

+ 1.32±0.62 

0.02328±0.02082 

(110)1:1 

(110)1:0 

2520.31 

2517.05 

-8.21±1.01 

-9.82±0.70 

-12.02±0.84 

+2.30±1.38 

+1.34±2.28 

-8.39±2.10 

0.04985±0.02419 


1793.00 

+3.90±0.33 

-0.83±0.41 

+7.38±0.28 

+11.27±0.54 

0.10104±0.01576 


PdaCu 


Face 

a 

Af 2 

Ag 

A23 

A& 

Ripple 


2438.15 

-13.86±3.92 

+0.48±0.41 

+8.50±0.41 


0.27452±0.08301 

1 1 mM 

1444.49 

-11.47±0.33 


+21.29±1.90 

-4.06±0.45 

0.21946±0.00637 

1 1 jtjrR 

2103.89 

-13.57±0.51 

-9.44±0.65 

+3.95±1.29 


0.05583±0.01565 

V f w 

2279.19 

-21.51±0.97 


+12.28±5.58 

-5.05±1.29 



1562.29 

-9.95±0.27 

-13.01±0.22 

-6.37±0.19 

-4.32±0.41 

0.06775±0.01091 


Cu3Pt 


I Face 

<T 

A$ 

A^ 

A23 

Ag 

Ripple 

B 1 2 1 SB 

2517.04 

-6.33±0.40 

+0.20±0.63 

+1.28±0.34 


0.12179±0.01928 

If W' 

2279.98 

+1.29±0.38 


+8.87±0.41 

+0.88±0.77 


1 f j 51 

2935.31 

-3.74±1.13 

+0.17±1.28 

+4.37±1.50 


0.05154±0.03176 

1 1 T 1 

2780.12 

-7.29±0.73 


+6.97±1.46 

-2.31 ±0.92 


RgftW 

1928.91 

+4.05±0.29 

+10.56±0.47 

+10.82±0.26 

+1.68±0.72 

0.14040±0.01644 


Table U Surface energies (in ergs/cm 2 ), planar relaxations (as percentages of the rigid 
interplanar spacing, see text) of the top two layers of several low-index faces of the LI 2 
ordered structure. A represents the relaxation between planes i and j of an atom of 
species X in layer i. The last column indicates the difference in position (in A) of an atom 
A and an atom B in the top layer. 
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Pt 3 Cu 


Face 

<r 

A £2 

*3 

A a 

^23 

A 23 

Ripple 

■ ! 

2002.83 

-6.12±0.35 

-7.70±0.31 

-3.00±0.25 


0.03136±0.01323 

III! 

1924.84 

2101.70 

-7.43±0.17 

-7.85±0.44 

-7.63rfc0.52 

-3.97±0.26 

-5.04±0.65 

-2.26±0.31 

0.00320±0.01364 

BB 

2583.47 

1146.30 

-4.63±0.10 

-8.61±0.15 

-3.27±0.15 

+1.51±0.12 

0.09149±0.00587 


AI 3 N 1 


Face 

a 


a£ 

A 23 

A® 

^23 

Ripple 

( 100 ) 1:1 

3435.00 

+5.31±2.02 

giifginggM 

+1.57±0.43 


0.21925±0.04706 

( 100 ) 1:0 

393.62 

-0.39±0.14 


+1.71±0.39 

+7.64±0.15 


( 110 ) 1:1 

2755.35 

-11.72±0.82 


-1.65±1.33 


0.00765±0.01948 

( 110 ) 1:0 

2260.07 

-1.42±0.86 


-10.55±1.07 

+15.64±0.94 



2410.59 

-8.04±0.19 


-6.14±0.16 

-1.34±0.23 

0.10936^0.00853 


N 13 AI 


Face 

a 

^12 

A® 

^12 

A 23 

A® 

^23 

Ripple 


2887.57 

-5.46±0.48 

+1.12±1.29 

+1.30±0.45 

H-0.00rfc0.00 

0.11819±0.03170 

1mm 

3156.69 

-0.54rfc0.49 


+5.94db0.56 

+0.66dtl.l4 


(110)1:1 

3169.63 

-7.17±1.24 

-3.00±2.26 

+0.29±2.34 


0.05310±0.04449 

luongl 

3885.99 

-7.98±0.79 


+3.08±1.92 

-1.40±2.16 


ImKal 

2575.54 

+5.08rfc0.39 

+13.05±1.18 

+11.07±0.38 

+2.34±1.13 

0.16545±0.03266 


AggAu 


| Face 

G 

Aft 

Ag 

A 23 

Al 

Ripple 

lj||l 

1537.59 

1488.54 

-20.70±3.46 

-5.35±0.49 

-6.07d:0.43 

-2.21dr0.41 

+13.85dr2.92 

-2.03±0.58 

0.30075rfc0.07993 

|i 

1661.68 

1620.23 

-6.71±0.75 

-11.38±1.04 

-5.16±0.72 

+0.81±1.37 

+2.68±5.09 

-3.43±1.15 

0.02262±0.02136 

(111)3:1 

1214.65 

-4.29dr0.35 

-4.84±0.28 

-1.74±0.27 

-1.89dr0.43 

0.01303±0.01484 


Au3Ag 


I Face 

G 

Aj2 

Ag 

A* 

^23 

Ag 

Ripple 


1578.82 

-4.29dr0.40 

-4.27rfc0.58 

-1.04±0.34 


0.00036dr0.02005 


1513.15 

-5.14±0.27 


-2.19dh0.44 

-1.21±0.48 

0.10575±0.00563 

(110)1:1 

1618.62 

-5.90rfc0.68 

-6.50±0.78 

-l.29drl.39 


0.00869±0.02118 


1657.45 

-5.61±0.42 


-0.67±1.47 

-0.94±1.80 


| (111)3:1 

1204.88 

-3.85±0.25 

-5.14dr0.43 

-1.50dr0.27 

-0.42dr0.52 

0.03046rfc0.01622 


Table 11 (Cont’d) 
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Au 3 Pd 


Face 

a 

^12 

^12 

^23 

7TB 

^23 

Ripple 


2301.77 

-1.5710.71 

-3.0610.72 

+0.0610.50 


0.0301910.02910 

(100)1:0 

1458.72 

-2.8810.34 


-0.5210.46 

+0.2910.53 


(110)1:1 

2138.84 

-4.3710.97 

-5.3811.07 

-0.1711.25 


0.0143710.02922 

(110)1:0 

1987.39 

-3.8510.66 


-0.5611.65 

+1.0111.85 


(111)3:1 

1838.26 

-3.4810.36 

-5.2610.50 

-1.6610.31 

-0.7310.69 

0.0418210.02015 


Pd 3 Au 


(100)1:1 

(100)1:0 

(110)1:1 

(110)1:0 

(111)3:1 

1584.61 

2304.90 

2104.72 

2318.14 

1357.94 

-26.59+2.83 

-3.22+0.72 

-4.79+1.16 

-10.22+1.26 

-0.60+0.56 

-1.02+0.60 

-2.88+1.17 

+7.46+0.44 

+1.26+0.57 

+14.53+3.37 

+0.93+1.67 

+7.92+5.15 

+8.77+0.44 

-1.12+0.71 

-2.94+1.10 

-0.19+0.52 

0.50770+0.06809 

0.06392+0.01425 

0.02681+0.03276 

0.18482+0.02292 

Ni 3 ] 

Pd 


2024.79 

-10.92+0.33 

-10.02+0.55 

-3.98+0.31 


0.01657+0.01635 


2513.61 

-9.24+0.36 


-7.63+0.42 

+0.40+0.44 


(110)1:1 

2477.88 

-9.20+0.68 

-13.05+0.56 

+3.19+0.98 


0.04994+0.01607 


2590.59 

-10.89+0.45 


+2.92+1.48 

-11.72+0.89 


(111)3:1 

1409.58 

+1.72+0.20 

+1.38+0.38 

+8.36+0.15 

-9.49+0.26 

0.00714+0.01237 


Pd 3 Ni 



a 

A 12 

AS 

A*^ 

^23 

/V s 

“23 

Ripple 

(100)1:1 

1966.28 

-19.69+3.37 

-16.82+0.26 

-8.30+0.21 


0.05543+0.06985 


1240.30 

-12.43+0.23 


+18.09+1.83 

-4.90+0.35 



1882.18 

-13.90+0.31 

-11.39+0.36 

+0.89+0.98 


0.03418+0.00922 


1809.03 

-22.77+0.67 


-13.57+0.38 

-6.29+1.06 


| (111)3:1 

617.93 

-10.77+0.13 

-14.22+0.09 

-7.07+0.07 

-4.99+0.24 

0.07678+0.00483 

Ni 3 ] 

Pt 

Face 

o 

&12 

AS 

A23 

AS 

Ripple 

(100)1:1 

2577.37 

-9.03+0.29 

-1.81+0.62 

-0.66+0.26 


0.13363+0.01679 

(100)1:0 

2759.94 

-0.46+0.29 


+7.99+0.30 

-0.30+0.71 


(110)1:1 

3355.64 

-4.79+0.77 

+3.49+1.32 

+8.72+1.08 


0.10835+0.02736 

(110)1:0 

3234.27 

-9.83+0.48 


+5.24+1.25 

-3.53+0.79 


(111)3:1 

2042.79 

+3.67+0.20 

+11.18+0.49 

+11.22+0.19 

+1.36+0.68 

0.16040+0.01477 


Table 11 (Cont’d) 


34 

















































Pt 3 Ni 


Face 

a 

A12 

Af, 

A 23 

Ag 

Ripple 


2885.91 

2100.63 

-1.61±0.69 

-3.55±0.23 

-5.80±0.38 

-0.70±0.34 

-1.47±0.41 

+1.63±0.32 

0.08166±0.02087 


2701.88 

2955.02 

-5.55±0.73 

-3.42±0.60 

-7.38±0.60 

-0.44±0.98 

-2.93±1.42 

+6.98±1.14 

0.02527±0.01832 

(111)3:1 

1823.75 

-3.18±0.15 

-8.05±0.18 

-2.37db0.16 

+2.84±0.17 

0.10961±0.00736 


Table ll(Cont’d) 
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Fig. 2 : Transformation seen by an atom ((g)) during the ideal process of alloy formation 
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Fig. 3: Breakup of transformation shown in fig. 2 into two independent processes: (a) strain 
energy and (b) chemical energy. The dots indicate atoms of species B while X indicates 
atoms of species A. <g> indicates the reference atom. 
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L1 0 (AB) LI, (AB) 


Fig. 4: Ordered structures of fee based alloys. 
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